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'Mass  resolved  excitation  vibronic  spectra  of  jet-cooled  7-azaindole  and  its  clusters  withal 

CUu  NH3,  H2O,  D2O,  CH3OR,  and  C2H5OH  axe  reported  and  analyzed  with  regard  to  molecular 

and  cluster  geometry  and  the  nature  of  the  first  excited  singlet  state.  Large  changes  in  the  various 

spectra  are  observed  upon  clustering  and  upon  deuteration  of  7-azaindole.  The  observed  vibronic 

spectra  of  both  7-azaindole  and  its  clusters  can  be  rationalized  with  two  general  assumptions: 

1.  the  hydrogen  attached  to  the  pyrrole  nitrogen  of  7-azaindole  is  out  of  the 

first  excited  singlet  state;  and  2.  the  observed  spectra  are  characterized 

mixing  not  completely  removed  by  the  clustering.  MOPAC  5  calculations  of  molecular  geometry 

:  i 

suggest  that  the  Sc*  state  is  nonplanar.  Additional  cluster  potential  energy  calculations  suggest  that 

the  formation  of  cyclic  hydrogen  bonded  clusters  is  not  likely  for  these  gas  phase  1:1  or  1:2  7; 

...  V ' 

azaindole/solvent  clusters:  the  major  gas  phase  solvent  clustering  probably  takes  place  at  the  ifc- 


molecular  plane  in  the 

-'V'  .  rr-liU 

by  strong  hit*  -Tlfif* 


system  of  the  7-azaindole  molecule.  Kj-  ■ 

1  ■  .  o  •  uv- 


V  "  r; 


]T 

■*  \  oA  C ; '  \  ..  • 
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I.  INTRODUCTION 

7-azaindole  (see  Figure  la,  7AZI)  is  suggested  to  have  two  veiy  interesting  and  important 
properties  in  solution:  1.  when  solvated  by  hydrogen  bonding  solvents  (such  as  water,  alcohols, 
ammonia,  or  even  itself),  it  is  thought  to  form  a  cyclic  structure  through  a  double  hydrogen  bond 
(Figure  lb)1’3;  and  2.  it  will  tautomerize  (Figure  lc)  in  protic  solvent  solutions  via  a  cyclic 
hydrogen  bonded  intermediate. 1>2'4’8  Such  behavior  is  suggestive  of  the  behavior  of  DNA  bases9, 
and  thus  7AZI  has  been  a  particularly  attractive  and  compelling  system  for  study. 

Spectra  of  isolated  7AZI,  and/or  7AZI  clustered  with  various  hydrogen  bonding  and  non¬ 
hydrogen  bonding  solvents  should  aid  in  the  elucidation  of  the  energetics  and  dynamics  involved  in 
both  the  cyclic  hydrogen  bonded  solvation  structures  and  in  the  tautomerization  process. 
Supersonic  jet  expansion  techniques  coupled  with  laser  spectroscopy  and  mass  resolution  (mass 
resolved  excitation  spectroscopy)  are,  of  course,  an  excellent  experimental  approach  to  the  study  of 
molecular  and  cluster  structure  and  energetics.  Coupled  with  state  of  the  art  molecular  orbital 
semiemperical  calculations  and  potential  energy  cluster  structure  calculations,  a  good  deal  can  be 
learned  about  bare  molecule  and  cluster  static  and  dynamic  properties.10’13 

Spectra  of  jet-cooled  7AZI,  (7AZI)2,  and  7AZI/water  clusters  have  been  previously 
obtained  with  a  view  toward  the  interpretation  of  the  tautomerization  process.14*16  The  cluster 
spectra  evidence  large  red  shifts  from  the  bare  molecule  spectrum  and  are  quite  complex:  in 
general,  the  bare  molecule  and  the  various  cluster  spectra  do  not  appear  to  coincide  with  regard  to 
either  vibronic  structure  or  intensity.  These  observations  were  attributed  to  double  hydrogen 
bonding  such  as  shown  in  Figure  1.  Nevertheless,  aside  from  this  explanation,  no  firm  evidence 
exists  for  such  structure  in  the  gas  phase.  In  fact,  evidence  from  other  theoretical  and  experimental 
studies  of  isolated  van  der  Waals  (vdW)  clusters10*13  strongly  suggests  that  the  major  interaction 
between  aromatic  (solute)  molecules  and  small  solvent  (i.e.,  H2O,  NH3,  CnH2n  +  2,  ROH,  etc.) 
molecules  is  through  the  Jt-electron  system  of  the  solute  and  the  major  electron  density  of  the 
solvent.  For  systems  with  a  substantial  hydrogen  bonding  interaction  in  addition  to  the  usual 
dispersion  interaction,  the  cluster  may  have  multiple  stable  structures:  a  planar  hydrogen  bonding 
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structure,  a  ring  centered  k -electron/electron  density  structure,  and  a  compromise  structure  in 
which  neither  form  of  interaction  dominates  the  geometry. 

In  addition  to  the  above  likely  complications  of  multiple  cluster  geometries  for  7AZI 
clusters  and  dimers,  the  isolated  molecule  may  evidence  its  own  complicated  spectral  behavior  due 
to  three  potential  causes:  1.  mixing  between  nearly  degenerate  'La  and  'Lb  tut*  excited  electronic 
states;  2.  mixing  between  the  lowest  kk*  state  (probably  'Lb  which  corresponds  to  the  'B2U  state 
of  benzene)  and  low  lying  nrc*  states;  and  3.  nonplanarity  of  the  7AZI  molecule  in  Si  -  in 
particular,  at  the  pyrrole  nitrogen.  The  'Lb  state  with  its  dipole  parallel  to  the  long  axis  of  the 
molecule  is  rather  insensitive  to  the  solvent  environment  while  the  'La  state  is  lowered  in  energy  by 
polar  environments.  In  the  case  of  7AZI,  both  transitions  are  polarized  more  or  less  along  the  long 
molecular  axis.17  The  'La  and  'Lb  states  are  well  separated  in  gas  phase  indole18  but  are 
apparently  close  together  in  7AZI.'*'9  For  most  isolated  indole-like  systems,  only  the  'Lb  state  is 

identified  at  low  energy,  mt*  -  7t7t*  interactions  have  not  been  reported  or  documented  for  7AZI 
but  have  been  suggested  for  purines20  and  other  systems21*22.  A  rotational  analysis  of  the  0° 

transition  for  Si  «-  So23  does  not  suggest  any  ntt*  -  kk*  mixing  but  almost  all  the  transition 
intensity  is  expected  to  be  of  a  kk*  character. 

The  work  reported  in  this  paper  addresses  the  above  three  complications  for  the  Si  state  of 
bare  and  clustered  7AZI.  Additionally,  structures  for  7AZI  clusters  with  argon,  methane,  water, 
alcohols,  and  ammonia  are  suggested  based  on  calculational  and  spectroscopic  data. 

II.  EXPERIMENTAL  PROCEDURES 

A  detailed  description  of  the  supersonic  jet/laser  apparatus  is  found  in  ref.  24.  The  doubled 
output  of  a  pulsed  Nd/YAG  laser  pumps  a  dye  laser  (R590,  R610,  KR620,  R640,  and  SR640  plus 
DCM)  whose  output  is  in  turn  doubled  into  the  ultraviolet  in  order  to  excite  the  samples  of  interest. 
The  spectral  range  covered  for  a  given  sample  often  includes  the  output  from  several  dyes  so  that 
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relative  intensities  must  be  calibrated  in  overlapping  dye  regions.  For  such  situations,  the 
intensities  are  noted  in  the  figure  captions. 

Mass  resolved  excitation  spectra  (time-of-flight-mass  spectroscopy-TOFMS)  are  obtained 
for  all  samples  discussed  in  this  report  as  the  7AZI  (Aldrich,  98%)  contains  impurities  that  strongly 
absorb  and  emit  light  in  the  spectral  regions  of  interest.  These  impurities  could  not  be  readily 
removed  by  recrystallization  or  vacuum  sublimation.  For  the  study  of  clusters,  2-color  TOFMS 
are  needed  to  obtain  clear  spectra  of  clusters  of  specific  stoichiometry. 

7AZI  is  placed  inside  the  head  of  a  pulsed  supersonic  nozzle  heated  to  ca.  70*C.  The 
expansion  is  generated  using  50-70  psi  of  helium  carrier  gas.  The  solvents  for  clustering  are  either 
mixed  with  the  carrier  gas  (gases)  or  placed  in  a  trap  before  the  nozzle. 

di-7AZI  (at  the  pyrrole  ring  nitrogen  -  Hjq  on  Ni,  see  Figure  1)  is  synthesized  by  mixing 
7AZI  (1%)  with  D2O  at  room  temperature,  extracting  the  di-7AZI  with  ether,  and  drying  the 
organic  component  in  vacuum.  NMR  spectra  taken  of  the  sample  confirm  that  only  the  most  acidic 
hydrogen  (H10  of  Figure  1)  is  replaced. 

III.  GEOMETRY  OPTIMIZATION  CALCULATIONS 

The  geometry  of  the  7AZI  molecule  is  calculated  for  the  So  and  Si  electronic  states  using 
the  Molecular  Orbital  Package  v.  5  (MOPAC  5)25.  MOPAC  5  uses  the  semiempirical  MNDO 
method  with  optimized  parametric  Hamiltonians  AMI26  and  PM327.  With  a  given  set  of  input  data 
for  bond  lengths  and  angles,  the  program  optimizes  molecular  geometry  and  point  charge 
distributions  over  the  atoms  to  obtain  the  minimum  heat  of  formation.  Additional  molecular 
constants  such  as  bond  characteristics,  dipole  moment,  and  normal  modes  of  vibration  can  also  be 
found.  These  calculations  are  carried  out  for  both  the  ground  and  first  excited  singlet  states  of 
7AZI. 

Two  methods  are  employed  to  calculate  the  minimum  energy  geometries  and  binding 
energies  for  7AZI/solvent  clusters:  the  above  MOPAC  5  program  and  a  potential  energy 


calculation  based  on  Lennard- Jones,  hydrogen  bonding,  and  coulomb  atom-atom  interactions.28 
These  latter  calculations  are  well  described  in  our  previous  publications.22 

Input  data  for  the  potential  energy  cluster  calculations  are  the  potential  parameters  for  the 
different  forms  of  interactions  and  different  atoms,  and  the  bond  angles,  bond  lengths  and  atomic 
charges  for  all  the  molecules.  These  molecular  data  come  from  MOP  AC  5  AMI  and  PM3 
Hamiltonian  calculations. 

IV.  RESULTS  AND  DISCUSSION 
A.  7-Azaindole  -  Spectroscopy 

The  1 -color  TOFMS  of  7AZI  is  presented  in  Figure  2.  The  intense  feature  at  34634  cm'1  is 
believed  to  be  the  Si  <-  So  transition  origin.  Low  energy  vibronic  features,  which  play  a  major 
role  in  our  overall  analysis  of  the  bare  molecule  and  cluster  spectra,  are  assigned  at  234,  280,  306, 
467,  and  483  cm'1.  Prominent  doublets  are  identified  in  this  spectrum  at  714, 718  cm-1  and  730, 
735  cm'1  and  an  intense  feature  is  located  at  934  cm*1 .  Upon  deuteration  of  7AZI  at  the  acidic  Hio 
position  (dj  -  7AZI)  two  changes  are  found  in  the  7AZI  spectrum:  1.  many  of  the  low  frequency 
modes  shift  to  lower  energy  by  as  much  as  3%  (Figure  2b);  and  2.  the  doublets  in  the  spectrum 
collapse  (Figure  3).  These  changes  are  tabulated  in  Table  I.  Both  of  the  above  observations  are 
quite  unusual  in  a  qualitative  and  quantitative  sense  for  a  "normal"  aromatic  Si  So  transition.29 
Additionally,  the  7AZI  spectrum  displays  a  large  number  of  vibronic  features  of  quite  weak 
intensity  which  are  not  very  obvious  from  the  figures.  This  is  another  striking  aspect  of  the  7AZI 
spectrum. 

Upon  comparing  the  7AZI  spectrum  with  that  of  indole18’30,  one  finds  few  common 
features.  The  vibronic  features  between  700  to  800  cm'1  are  similar  for  both  molecules  but  the 
doublet  structure  here  and  throughout  the  spectrum  (for  weaker  features)  is  unique  to  7AZI.  The 
low  energy  peaks  at  234,  280,  306, 467,  and  483  cm'1  are  also  unique  to  7AZI.  These  differences 
between  the  two  spectra  are  surprising  considering  the  Si  <-  So  transition  is  supposed  to  be  to  the 


same  nn*  *Lb  Si  state  in  both  instances.  Moreover,  indole  has  only  two  weak  vibronic  features 
below  600  cm*1 18*20,  while  7AZI  has  more  than  13  relatively  intense  low  energy  features. 

The  unusually  large  number  of  low  lying  vibronic  features  for  7AZI  and  the  difference 
between  the  7A23,  di-7AZI  and  indole  spectra  can  be  accounted  for  by  two  molecular  properties  of 
7AZI:  1.  the  molecule  is  nonplanar  in  Si,  particularly  at  the  Ni-Hiopyrrolic  site;  and  2.  strong 
me*  -  TtJt*  vibronic  mixing  occurs  in  the  region  of  the  Si  <—  So  transition. 

The  nonplanarity  of  7AZI  in  the  Si  state  (in  particular  at  Hio)  would  account  for  the  large 
deuterium  substitution  (shifts  and  collapse  of  doublets)  effect  on  the  spectrum  and  might  also  lead 
to  an  increased  number  of  allowed  vibronic  transitions.  Even  this  nonplanarity,  however,  should 
not  give  rise  directly  to  so  many  transitions. 

We  suggest  that  a  strong  mi*  -  nil*  mixing  would  account  for  the  apparent  density  of 
vibronic  transitions  in  the  Si  <-  So  manifold.  Similar  suggestions  have  been  made  for 
isoquinoline21  ^  1,  quinoline21, 2-hydroxypyridine22  and  other  molecules.  The  intensity  of  these 
vibronic  transitions  comes  from  the  nn*  state  but  the  peak  positions  are  associated  with  mr*  -  nn* 
vibronic  resonances  and  couplings.  As  in  the  case  of  isoquinoline,  we  arc  unable  to  detect  direct 
absorption  (through  either  TOFMS  or  FE  experiments)  from  So  to  the  mt*  excited  state. 

In  order  to  elucidate  this  problem  further,  two  approaches  are  useful:  1.  calculation  of  the 
So  and  Si  geometry;  and  2.  clustering  of  7AZI  with  various  solvents  to  alter  the  nn*  -  nn* 
coupling  and  perhaps  to  affect  the  7AZI  planarity. 

B .  7-Azaindole  -  Calculations 

The  geometry  and  atomic  point  charge  distribution  for  7AZI  in  both  the  ground  and  first 
excited  singlet  states  are  calculated  using  MOPAC  5  AMI  and  PM3  Hamiltonians.  These  data  are 
summarized  in  Tables  II  and  HI.  The  optimum  geometries  calculated  with  either  Hamiltonian  are 
the  same  for  the  ground  state,  although  significant  differences  are  found  for  the  atomic  partial 
charges  at  various  atoms  (e.g.,  C9  and  Ni).  For  the  first  excited  state,  the  two  Hamiltonians  give 
different  geometries:  the  AMI  Hamiltonian  yields  a  nonplanar  molecule  (C9  and  Hio)  while  the 
PM  3  Hamiltonian  yields  a  planar  molecule.  These  calculations,  however,  arc  neither  very  accurate 


nor  reliable  as  the  convergence  criteria  for  optimum  Si  geometry  for  both  these  calculations  is  poor 
(»  7  for  Si  as  opposed  to  0.01  for  So).  Moreover,  different  optimum  geometries  obtain  for 
different  initial  input  geometries.  The  general  trend  of  results  for  the  AMI  calculations  seems 
sound,  however,  7AZI  apparently  has  a  propensity  for  nonplanarity  of  the  pyrrolic  ring  in  Si, 
especially  at  C9  and  H10.  The  AMI  Hamiltonian  can  generate  a  nonplanar  7AZI  for  the  first 
excited  singlet  state.  H10  is  typically  out  of  plane  by  18-19*  (see  Figure  4). 

We  can  conclude  from  both  the  experimental  and  theoretical  results  that  some  of  the  spectral 
crowding  and  doubling  of  various  vibronic  features  in  the  7AZI  spectrum  may  well  be  due  to  the 
loss  of  planarity  in  Si.  The  potential  well  for  H10  motion  with  respect  to  the  plane  of  the  heavy 
atom  ring  is  probably  a  double  well  with  a  large  barrier  to  the  zero  degree,  planar  conformation. 

C.  7-Azaindole(Ar)i  and  (CH4)i  clusters 

The  general  reason  to  look  at  cluster  spectra  at  this  point  is  to  resolve  the  two  issues  posed 
above:  1.  the  (Si)  nonplanarity  of  7 AZI  and  its  ensuing  double  well  potential  surface;  and  2.  the 

importance  of  ntt*  -  jnc*  mixing  for  this  molecule. 

Spectra  for  7AZI(Ar)i  and  (CHOi  are  presented  in  Figure  5.  The  cluster  shifts  for  the  0q 

transition  are  -18  cm'1  for  7AZI(Ar)i  and  -46  cm*1  for  7AZI(CH4)i.  These  are  quite  typical  of 
cluster  shifts  for  other  solute  systems.  For  example  for  toluene,  pyrazine,  pyrimidine,  benzene, 
and  indole  the  origin  shift  for  1:1  methane  clusters  are  -4312,  -3313,  -5732,  -41 H,  and  -3633  cm*1, 
respectively  and  the  shifts  for  benzene  and  indole  with  argon  are  -21 34  and  -2635  cm*1, 
respectively. 

Calculations  of  cluster  geometry  are  presented  in  Figure  6.  Again  very  typical  behavior  is 
observed  for  the  single  minimum  energy  configurations  of  these  clusters. 

By  contrast,  comparison  of  the  vibrational  structure  of  these  clusters  with  each  other  and 
bare  7AZI  shows  substantial  changes  in  the  vibronic  features  brought  about  by  clustering.  This  is 
strong  evidence  that  n n*  -  kk*  mixing  is  distorting  the  observed  cluster  and  bare  molecule  vibronic 
structure.  Thus  the  apparent  density  of  vibronic  features  in  the  7 AZI  bare  molecule  spectrum  may 
well  be  due  to  vibronic  coupling  between  a  low  lying  ntt*  state  and  the  "observed"  kk*  state,  "Si". 


The  observed  vibronic  transitions  in  this  region  can  reflect  this  interaction  and  not  existing 
vibrational  structure  in  the  7AZI  excited  xx*  electronic  state. 

At  higher  vibrational  levels  in  Sj  the  transition  intensity  of  both  clusters  begins  to  decrease 
due  to  vibrational  predissociation  of  the  clusters.36  The  intense  feature  at  +934  cm*1  above  the  Oq 

transition  of  the  bare  molecule  is  no  longer  observed  either  due  to  vibrational  predissociation  or 
removal  of  the  nx*  -  xx*  interaction  upon  clustering.  In  either  event,  this  demonstrates  that  the 
934  cm*1  feature  is  a  vibronic  one  and  not  a  new  electronic  origin  (e.g.,  'La ).  Thus  the  observed 
spectra  are  suggested  to  be  due  to  a  transition  involving  the  'Lb  (xx*).  The  excited  'Lb  (rat*)  state 
probably  interacts  with  a  lower  energy  singlet  nx*  state  to  produce  the  observed  complex  high 
density  Si  <-  So  transition  in  the  34,600  to  35,600  cm*'  region. 

D.  7-Azaindole/H20  Clusters 

The  spectra  of  7AZI(H20)ii2,3  are  presented  in  Figure  7.  The  Oq  transition  for 
7AZI(H20)i  is  located  at  33,346  cm*',  red  shifted  by  1288  cm*'  from  the  bare  molecule  Oq 

transition.  The  7AZI(H20)2,3  spectra  are  further  red  shifted  by  1997  cm*1  and  2075  cm*1, 
respectively.  These  red  shifts  are  enormous  considering  those  found  for  benzene32  (+85  cm*'), 
indole37  (-135  cm*'),  phenol38  (-  356  cm*'),  fluorene39  (+  57  cm*'  ),  and  isoquinoline2' 
(0.6  cm-')  clusters  with  (H20)].  The  unusually  large  red  shift  for  7AZI/water  clusters  suggests 
that  the  nx*  -  7t7t*  excited  state  interaction  has  been  (at  least  partially)  removed  by  the  interaction  of 
7AZI  with  water.  Hydrogen  bonding  solvents  are  known  to  raise  the  energy  of  nx*  states  with 
respect  to  xx*  states.40  The  suggested  mechanism  for  this  increase  in  the  nx*  excitation  energy  is 
the  lowering  of  the  ground  state  n  orbital  energy  due  to  hydrogen  bond  formation  between  water 
and  the  lone  pair  electrons  on  the  ring  nitrogen  (see  Figure  lb).  If  the  nx*  state  were  raised  in 
energy  above  the  xx*  state,  additional  lowering  in  energy  of  the  xx*  state  would  occur  through 
interaction  between  these  two  zero  order  descriptions. 

Not  only  is  the  origin  shift  of  7AZI/water  clusters  with  respect  to  7AZI  suggestive  of 
excited  electronic  state  mixing  but  the  observed  cluster  vibronic  structure  is  as  well.  Note  that  the 


low  energy  modes  of  /AZI  Si  have  changed  nearly  33%  in  7AZI(H20)i:  this  is  almost  unique  for 
vdW  clusters.  Apparent  overtone  structure  is  observed  for  both  systems.  The  strong  7AZI  feature 
at  934  cm*1  is  missing  in  7AZI(H20)n  spectra;  vibrational  predissociation  is  unlikely  to  be  the 
reason  for  this  absent  feature  based  on  binding  energy  considerations. 

Given  the  large  cluster  shifts,  the  change  in  Si  vibrational  structure  upon  clustering  with 
water,  and  the  absence  of  any  other  7AZl/water  features  in  the  33,000  to  36,000  cm*1,  we  believe 
that  the  excited  state  observed  for  7AZI/water  is  'Lb  with  a  significant  removal  of  mt*  state 
interference.  The  'La  state  is  probably  not  observed  in  this  region. 

Cluster  structure  has  been  explored  through  a  number  of  different  calculations.  Both 
MOPAC  5  and  potential  energy  calculations  are  performed  on  the  7AZI(H20)i  clusters.  Table  IV 
and  Figures  8  and  9  summarize  the  results  of  such  studies.  The  MOPAC  5  calculation  gives 
distorted  double  hydrogen  bonded  structures  for  the  PM3  Hamiltonian  (Figure  8a)  but  single 
hydrogen  bonded  structures  for  the  AMI  Hamiltonian  (Figure  8b).  In  neither  case  are  the  water  or 
7AZI  structures  and  charges  much  different  from  those  of  the  isolated  individual  molecules.  Note 
that  the  MOPAC  5  calculation  [tested  in  this  study  for  the  pyrazine41  and  pyrimidine4!  dimers  and 
pyrazine(NH3)i32]  tends  to  underestimate  the  vdW  dispersion  interaction  between  the  aromatic  7t- 
system  and  the  solvent  heavy  atom(s).  In  general  these  structures  do  not  appear  to  be  very  reliable. 

The  potential  energy  calculation  (Figure  9)  generates  different  structures  for  7AZI(H20)i 
depending  on  the  Hamiltonian  employed  for  the  isolated  molecule  structure  calculations.  The  PM3 
molecules  generate  a  single  structure  while  the  AMI  molecules  yield  three  different  structures,  one 
of  which  is  the  PM3  structure.  Nevertheless,  all  of  these  structures  are  the  basic  compromise 
structure  falling  between  those  found  for  only  hydrogen  bonding  interactions  and  those  found  for 
only  dispersion  interactions.  If  a  double  hydrogen  bonded  structure  (Figure  8a)  is  fixed  for  the 
potential  energy  calculation  as  an  initial  geometry,  a  binding  energy  of  ca.  560  cm*1  obtains;  this 
structure  is,  however,  not  at  a  minimum  on  the  potential  surface.  Such  minimum  energy  structures 
are  presented  in  Figure  9.  This  of  course  makes  a  good  deal  of  "intuitive  sense".  We  believe  that 
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these  geometries  are  more  reliable  than  either  the  PM3  or  AMI  MOPAC  5  cluster/dimer  geometries 
(Figure  8).  , 

Potential  energy  calculations  for  higher  order  clusters,  based  on  PM3  molecules,  generate 
two  significant  results:  a  plethora  of  structures,  and,  if  the  ring  center  positions  are  occupied,  a 
double  hydrogen  bonded  cyclic  structure.  The  suggestion  here  is  that  the  largest  interaction 
between  7AZI  and  water  is  at  the  rc- system  of  7AZI,  but  that  once  these  positions  are  occupied  the 
hydrogen  bonding  structures  (cyclic  and  otherwise)  can  be  realized. 

E.  7-Azaindole/D20  Clusters 

If  D2O  is  substituted  for  H2O  in  the  expansion  system  mixture,  mass  resolved  excitation 
spectra  are  observed  for7A23  (118  amu),  (118  +  1)  amu,  (118  +  18)  amu,  (118  +  19)  amu,  (118 
+  20)  amu,  and  (118  +  21)  amu,  etc.  The  spectra  at  mass  channels  (118  +  1)  and  (118  +  18)  amu 
are  identical  to  those  of  di  -  7AZI  and  7AZI(H20)i,  respectively.  Spectra  taken  at  mass  channel 
(118  +  18),  (118  +  19),  (118  +  20)  and  (118  +  21)  amu  are  presented  in  Figure  10.  The  spectra 
certainly  suggest  that  a  rapid  exchange  between  the  acidic  H10  of  7AZI  and  the  protons  of  water 
takes  place  on  the  So  potential  jrface.  Thus  the  spectrum  of  Figure  10b  is  composed  of  dj  - 
7AZI(H20)i  and  7AZI(HOD)i  spectra  and  the  spectrum  of  Figure  10c  is  composed  of  dj  • 
7AZI(HOD)i  and  7AZI(D20)i  spectra.  The  spectrum  taken  at  mass  channel  (118  +  21)  amu 
(Figure  lOd)  is  then  due  to  di  -  7AZI(E>20)i.  The  ratio  of  peak  heights  in  these  spectra  is  roughly 
that  expected  for  the  statistical  distribution  of  the  various  isotopic  cluster  species. 

Two  features  of  Figure  10  are  worth  special  emphasis.  First,  both  spectra  at  (1 18  +  19) 
amu  and  (1 18  +  20)  amu  consist  of  three  origin  transition:  three  origins,  of  course,  implies  the 
existence  of  three  distinct  sites  for  the  deuterium  substitution.  The  three  sites  are  most  likely  H10 
of  7AZI  (as  found  specifically  for  this  particular  deuteration)  and  the  two  distinct  proton  sites  on 
the  water  molecule  in  the  cluster.  This  experimental  observation  is  not  consistent  with  the  AMI 
MOPAC  5  geometry  depicted  in  Figure  8b.  Second,  the  spectrum  of  di  -  7AZI(D20)i  is  much 
simpler  (fewer  transitions)  than  that  of  7AZI(H20)i  (compare  Figures  10a  and  lOd).  The  nn*  - 
7t7t*  interaction  suggested  above  for  7AZI  may  thus  not  be  totally  removed  even  in  the  presence  of 
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water  solvation.  This  is  perhaps  further  evidence  against  a  strong  hydrogen  bonded  cyclic 
structure  for  7AZI(H20)i  as  generated  by  the  PM3  MOP  AC  5  cluster  calculation  (Figure  8a). 

Figure  1 1  shows  the  spectrum  taken  at  mass  channel  (1 18+37).  The  mass  corresponds  to 
7AZI  plus  two  waters  with  one  of  the  water  hydrogens  or  Hio  is  dueterated.  The  five  possible 
conformers  due  to  the  different  dueterium  substitution  sites  appear  in  three  distinctive  bands  of 
electronic  origins.  The  band  B  apparently  represents  a  single  peak  and  can  be  tentatively  assigned 
to  di-7AZI(H20)2,  while  each  of  bands  A  and  C  is  apparendy  composed  of  two  very  close  peaks 
and  can  tentatively  be  assigned  to  7AZI(H20)(DOH).  Attempts  to  identify  the  structure  of  7AZI 
plus  two  waters  from  Figure  1 1  meet  with  some  difficulty;  however,  we  believe  that ,  if  the  double 
hydrogen  bonding  structure  in  Figure  8a  were  realistic,  four  distinctive  bands  would  be  observed 
at  mass  channel  (118+37)  for  the  cluster  origin. 

F.  7-Azaindole/Alcohol  Clusters 

Figure  12  shows  the  2-color  TOFMS  of  7AZI(CH30H)i  and  7AZI(C2H50H)i. 
7AZI(CH30H)i  displays  a  single  origin  while  7AZI(C2HsOH)i  shows  multiple  origins.  The 
spectra  do  not  display  common  vibronic  features,  again  suggesting  that  the  mt*  -  tcjc*  interaction 
has  not  been  removed  in  these  clusters. 

The  potential  energy  calculations  for  7AZI(CH30H)j  are  presented  in  Figure  13.  Two 
structures  are  found  for  this  system.  The  hydrogen  bonding  interaction  weakens  for  ROH 
compared  to  HOH  and  the  electron  density  of  R  is  far  greater  than  that  of  H;  the  dispersion 
interaction  between  the  solvent  electron  density  and  the  ring  Tt-system  becomes  more  significant  in 
these  instances.  The  ethanol  cluster  has  many  different  calculated  conformations  due  mostly  to 
orientations  of  the  ethyl  group:  basically  the  cluster  geometry  for  7AZI(CH30H)i  and  (C2HsOH)i 
clusters  are  similar. 

G .  7-AzaindoIe/Ammonia  Clusters 

7AZI  clustered  with  ammonia  has  a  spectrum  different  from  both  the  bare  molecule  7AZI 

and  the  other  clusters,  as  presented  in  Figure  14.  Note  especially  the  low  lying  vdW  modes 
following  within  50  cm*1  of  the  0q  transition  and  the  low  lying  molecular  vibronic  features  at 
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ca.  200  and  400  cm*1.  The  7AZI(NH3)i  spectrum  looks  much  simpler  (i.e.,  has  fewer 
transitions)  than  the  7AZI(H20)i  spectrum. 

According  to  the  potential  energy  calculation,  the  NH3  molecule  coordinates  to  the  aromatic 
Jt  -  system  for  the  minimum  energy  configuration  (See  Figure  15).  The  7AZI(NH3)2  cluster  has 
the  typical  symmetric  and  asymmetric  structures  as  depicted  in  Figure  15.  The  observed  spectrum 
for  7AZI(NH3)2  is  thought  to  arise  from  the  symmetric  structure  because  of  the  cluster  shift,  the 
large  number  (five)  of  different  asymmetric  structures  generated  in  the  calculation,  and  the  fact  that 
80%  of  the  calculated  minimum  energy  structures  are  of  the  form  displayed  in  Figure  15b. 

H .  van  der  Waals  Modes  for  the  7AZI  Ousters 

Figure  16  summarizes  the  vdW  mode  spectra  for  the  various  clusters.  Note  that  in  general 

the  intensity  of  these  features  is  weak  and  the  Franck-Condon  envelop  for  the  transition  has  its 

maximum  at  the  Av  =  0  position.  Similar  results  are  obtained  for  most  vdW  clusters.31  None  of 
the  features  above  0q+  100  cm*1  in  the  various  clusters  can  be  attributed  to  vdW  modes. 

I.  Ouster  Shifts 

Ouster  red  shifts  are  large  and  probably  more  indicative  of  the  njt*-7i7i*  residual  mixing  in 
the  cluster  than  they  are  of  solute/solvent  interaction  energy.  For  example,  the  shifts  presented  in 
Table  V  do  not  reflect  hydrogen  bonding  strengths.  Moreover,  the  addition  of  a  second  solvent 
molecule  also  generates  a  large  red  shift  in  the  disolvate  cluster.  Both  observations  seem  to 
suggest  that  a  planar  double  hydrogen  bonded  structure  is  not  the  lowest  energy  structure  for  these 
solute/solvent  systems  in  the  gas  phase. 

V.  SUMMARY  AND  CONCLUSIONS 

Results  of  these  studies  for  7AZI  and  its  clusters  with  argon,  methane,  water,  alcohols,  and 
ammonia  can  be  summarized  as  follows:  1.  the  spectrum  of  the  bare  molecule  is  quite  congested 
with  many  low  energy  vibronic  features  and  apparent  doublets  appearing  in  the  mass  resolved 
excitation  spectrum;  2.  the  cluster  spectra  are  also  of  the  same  nature  but  the  vibronic  structure  is 
quite  different  for  each  particular  cluster,  3.  the  vibronic  doublets  disappear  upon  deuteration  of 
7AZI  at  H10;  4.  calculations  suggest  that  7AZI  may  be  nonplanar  in  Sj;  5.  potential  energy 
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calculations  suggest  that  cyclic  hydrogen  bonding  between  water,  alcohols  and  ammonia  and  7AZI 
is  not  a  stable  low  energy  structure  for  the  first  two  solvent  molecules  in  the  cluster,  and  6.  all 

j 

cluster  spectra  are  red  shifted  from  the  bare  molecule  origin,  especially  for  polar  solvents  as 
indicated  in  Table  V. 

The  conclusions  we  draw  from  these  results  are  as  follows:  1.  strong  rnt*  -  rot*  vibronic 
mixing  occurs  between  the  *Lb  rot*  excited  state  and  a  lower  lying  nn*  state;  2.  this  vibronic 
mixing  is  only  partially  removed  by  clustering;  3.  a  cyclic  hydrogen  bonded  structure  is  not 
achieved  for  the  7AZI(H20)i,2  7AZI(ROH)it2,  or  7AZI(NH3)i,2  clusters;  and  4.  the  7AZI 
molecule  is  nonplanar  (at  Hjo  -  Ni  on  the  pyrrole  ring)  in  the  first  excited  singlet  state. 

The  large  spectral  S  i  <—  So  red  shifts  for  water,  alcohol,  and  ammonia  7AZI  clusters  do  not 
necessarily  imply  a  cyclic  hydrogen  bonded  cluster  structure  because  the  (supposed)  stability  of 
these  structures  would  not  correlate  with  size  of  the  red  shift. 

Finally,  these  results  suggest  that  the  double  hydrogen  bonding  structures  assumed  to  be 
intermediates  in  the  7AZI  tautomerization  process  in  condensed  phases  may  only  be  present  in 
rather  high  solvent  density  systems  for  which  solvent  molecules  occupy  the  ring  centered  solvation 
sites  of  7AZI.  Under  such  conditions,  solvent  molecules  may  occupy  the  lower  binding  energy 
hydrogen  bonded  sites  at  the  N7  and  N1-H10  positions  (See  Figure  la). 
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Table  I 

Deuterium  substitution  (Hio,  see  Figure  1)  effect  on  vibrational  energies  of  7AZI.  (Not  all  weak  peaks 
after  700  cm*l  are  tabulated.) 


Vibrational  energies  (cm**) 

Intensity 

7AZ3 

dl-7AZI 

0(34634) 

0(34637) 

vs,  0q  origin 

234 

230 

m 

280 

272 

s 

306 

301 

m 

372 

366 

w 

404,  409 

397,400 

w 

457 

453 

w 

467 

458 

m 

483 

475 

s 

527 

523 

w 

547 

545 

w 

574,  578 

567,  573 

w 

672 

652 

w 

682 

688 

w 

696 

695 

w 

706 

700 

w 

714,718 

708,  709 

m 

730,  735 

725,  726 

s 

747 

735, 736 

m 

757 

742,  743 

s 

772 

750,752 

w 

796 

m 

857,  866 

851,  859 

w 

934 

934 

s 

16 


Table  li 

Optimum  geometry  and  charge  distribution  for  7AZI  in  the  ground  state.  Results  from  MOP  ACS 
calculations  using  the  PM3  and  AMI  Hamiltonians. 


^tom«* 

I 

Bond  length 

(A) 

PM3  AMI 

Bond  Angle 

(degrees) 

PM3  AMI 

Twist  Angle 

(degrees) 

PM3  AMI 

J 

K 

Charge 

L*>  (e) 

PM3  AMI 

1  N 

0.3155 

-0.2033 

2C 

1.410 

1.397 

1 

-0.2419 

-0.0762 

3  C 

1.381 

1.392 

109.27 

110.59 

2 

1 

-0.1466 

-0.1921 

4  C 

1.403 

1.396 

134.72 

135.68 

-179.93 

180.00 

9 

3 

2 

-0.0035 

-0.0420 

5  C 

1.378 

1.386 

117.62 

117.99 

179.98 

180.00 

4 

9 

3 

-0.1871 

-0.2116 

6  C 

1.417 

1.419 

120.96 

119.86 

0.01 

0.03 

5 

4 

9 

-0.0454 

-0.0530 

7  N 

1.337 

1.336 

123.17 

125.34 

-0.02 

-0.02 

6 

5 

4 

-0.0687 

-0.1454 

8  C 

1.402 

1.402 

107.95 

108.14 

0.15 

0.09 

1 

2 

3 

-0.1681 

0.0281 

9  C 

1.424 

1.459 

107.87 

107.59 

0.12 

-0.07 

8 

1 

2 

-0.1300 

-0.1210 

10  H 

0.987 

0.985 

125.66 

126.46 

179.55 

179.87 

1 

2 

3 

0.0768 

0.2592 

11  H 

1.090 

1.092 

128.40 

129.33 

179.98 

179.98 

2 

3 

9 

0.1400 

0.1676 

12  H 

1.087 

1.085 

126.77 

126.79 

179.91 

179.96 

3 

2 

1 

0.1285 

0.1587 

13  H 

1.094 

1.099 

120.50 

120.58 

180.01 

179.99 

4 

9 

8 

0.1061 

0.1386 

14  H 

1.094 

1.097 

120.43 

120.95 

180.00 

180.02 

5 

4 

9 

0.1118 

0.1407 

15  H 

1.096 

1.106 

121.34 

119.75 

179.99 

179.99 

6 

5 

4 

0.1125 

0.1568 

Final 

neat  of  formation:  48.15  kCal/mol  for  P 
66.97  kCal/mol  for  A 

V13 

Ml 

Footnotes: 

a:  Numbering  of  atoms  are  as  specified  in  Figure  la. 

b:  The  bond  length  is  defined  between  I  and  J. 

The  bond  angle  is  defined  between  lines  (U)  and  (JK). 

The  twisting  angle  is  defined  between  plane(JKL)  and  line  (LI). 
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Table  III 

Optimum  geometry  and  charge  distribution  for  7AZI  in  the  first  excited  singlet  state.  Results  from 
MOPAC5  calculations  using  the  PM3  and  AMI  Hamiltonians. 


Atom 

I 

Bond  length 

(A) 

PM3  AMI 

Bond  Angle 
(degrees) 

PM3  AMI 

Twist  Angle 
(degrees) 

PM3  AMI 

J 

K 

L 

Charge 

(e) 

PM3  AMI 

IN 

0.4769 

-0.1686 

2C 

1.364 

1.362 

1 

-0.4079 

-0.1737 

3C 

1.436 

1.459 

108.34 

110.04 

2 

1 

-0.0031 

-0.0743 

4C 

1.385 

1.404 

134.98 

134.96 

-179.80 

179.85 

9 

3 

2 

-0.1176 

-0.1329 

5C 

1.418 

1.404 

117.18 

116.89 

179.82 

-179.88 

4 

9 

3 

-0.1255 

-0.1619 

6C 

1.377 

1.395 

121.63 

120.82 

0.08 

1.22 

5 

4 

9 

-0.0854 

-0.0926 

7N 

1.395 

1.365 

122.88 

125.20 

0.01 

-0.42 

6 

5 

4 

-0.1268 

-0.1743 

8C 

1.422 

1.446 

109.76 

108.23 

-0.18 

-3.57 

1 

2 

3 

-0.1590 

0.0887 

9C 

1.422 

1.454 

107.31 

107.23 

0.20 

4.20 

8 

1 

2 

-0.1307 

-0.1330 

10H 

0.986 

0.993 

125.28 

126.47 

179.86 

-161.22 

1 

2 

3 

0.0648 

0.2510 

11H 

1.092 

1.092 

126.66 

127.26 

179.95 

-176.44 

2 

3 

9 

0.1585 

0.1760 

12H 

1.088 

1.085 

125.65 

125.26 

179.95 

-177.45 

3 

2 

1 

0.1252 

0.1630 

13H 

1.092 

1.095 

121.73 

121.37 

179.93 

177.86 

4 

9 

8 

0.1111 

0.1393 

14H 

1.094 

1.099 

118.59 

119.66 

180.03 

-178.98 

5 

4 

9 

0.1047 

0.1267 

15H 

1.094 

1.105 

123.41 

121.08 

180.02 

-178.85 

6 

5 

4 

0.1147 

0.1562 

Final  Heat  of  Formation: 

145.81  kCal/mol  for  PM3 
170.05  kCaVmol  for  AMI 

Table  IU 


Optimum  geometry  and  charge  distribution  for  7AZI(H20)i  in  the  ground  state.  Results  from  MOPAC5 
calculations  using  the  PM3  and  AMI  Hamiltonians. 


\tom 

I 

Bond  length 

(A) 

Bond  Angle 
(degrees) 

Twist  Angle 
(degrees) 

J 

K 

L 

Charge 

(e) 

1  N 

0.3152  -0.2024 

2  C 

1.409 

1.394 

1 

-0.2368  -0.0710 

3  C 

1.381 

1.392 

109.35 

110.57 

2 

1 

-0.1505  -0.1988 

4  C 

1.398 

1.396 

134.82 

135.89 

■179.70 

179.94 

9 

3 

2 

0.0014  -0.0441 

5  C 

1.379 

1.387 

117.63 

118.01 

■179.86 

-179.89 

4 

9 

3 

0.1900  -0.2133 

6  C 

1.416 

1.418 

120.96 

119.89 

0.09 

0.03 

5 

4 

9 

-0.0428  -0.0517 

7  N 

1.338 

1.336 

123.22 

125.30 

-0.05 

0.00 

6 

5 

4 

-0.0882  -0.1728 

8  C 

1.402 

1.402 

107.88 

108.04 

0.63 

0.01 

1 

2 

3 

-0.1692  0.0232 

9  C 

1.427 

1.460 

108.05 

107.71 

-0.30 

0.02 

8 

1 

2 

-0.1299  -0.1189 

10H 

0.989 

0.993 

125.90 

126.89 

179.34 

-179.90 

1 

2 

3 

0.0930  0.2809 

11  H 

1.090 

1.092 

128.40 

129.26 

179.98 

180.00 

2 

3 

9 

0.1410  0.1692 

12  H 

1.087 

1.085 

126.80 

126.81 

179.80 

180.03 

3 

2 

1 

0.1288  0.1586 

13  H 

1.094 

1.099 

120.52 

120.59 

179.97 

180.01 

4 

9 

8 

0.1064  0.1391 

14  H 

1.094 

1.097 

120.45 

120.95 

180.02 

179.98 

5 

4 

9 

0.1124  0.1409 

15  H 

1.096 

1.106 

121.28 

119.68 

179.89 

180.00 

6 

5 

4 

0.1105  0.1540 

160 

3.352 

3.083 

150.22 

156.81 

•174.78 

179.76 

7 

6 

5 

-0.3851  -0.4082 

17  H 

0.953 

0.962 

33.93 

66.66 

-44.12 

-61.76 

16 

7 

6 

0.2002  0.2075 

18  H 

0.951 

0.962 

107.61 

103.36 

-86.61 

-49.24 

16 

17  7 

0.1835  0.2075 

1  Final  heat  of  formation:  -7.71  kCal/mol  for  PM3 

1.04  kCal/mol  for  AMI 

Table  U 


Solvent  shifts  for  0®  electronic  transition  energy  of  1AZL 


System 

Cluster  shift 
(cm'1) 

Cluster  binding  energy3, 

(cm'1) 

7AZI 

0 

__ 

7  AZI(  Ar)  l 

-18 

-345 

7AZI(CH4)1 

-46 

-645 

7AZI(H20)i 

-1288 

-819 

7AZI(H20)2 

-1997 

-1630,  -2760, -2900 

7AZI(H20)3 

-2075 

-3500  -  -4500 

7AZI(CH30H)i 

-1425 

-1000,-1100 

7AZI(CH30H)2 

-2376*> 

7AZI(C2H50H)i 

-1454,-1420,-1430 

-1000- -1200 

7AZI(NH3)i 

-1217 

-930 

7AZI(NH3)2 

-1694 

-1850,-1950 

(7AZI)2C 

-2342 

-2310 

a:  based  on  PM3  Hamitonian  to  optimize  7AZI  charges  and  geometry, 
b:  reference  16. 

C:  not  shown  in  this  paper. 


20 


Figure 

Figure  1: 

Figure  2: 

Figure  3: 
Figure  4: 
Figure  5: 
Figure  6: 

Figure  7: 
Figure  8: 
Figure  9: 

Figure  10: 
Figure  11: 
Figure  12: 
Figure  13: 
Figure  14: 


Captions 

The  structures  of  (a)  7AZI  (b)  7AZI/water  cluster  with  a  double  hydrogen  bonding 
structure  and  (c)  (7AZI)2  and  its  possible  tautomerization.  The  numbering  on  atoms  in 
a)  is  used  in  describing  atoms  in  text  and  in  Tables  II,  El,  IV. 

1-color  TOFMS  of  a)  7AZI  and  b)  di*7AZI  (Hio)  for  the  region  0°  +1000  cm-*.  The 

peak  at  34634  cm*l  is  assigned  as  the  7AZI  origin.  The  origin  features  is  2  or  3  times 
larger  than  shown. 

1- color  TOFMS  of  a)  7AZI  and  b)  di-7AZI  with  an  expanded  energy  scale.  Note  the 
collapse  of  7AZI  doublet  structure  in  this  region  upon  Hjo  dueteration. 

MOP  AC  5  result  for  the  first  excited  singlet  state  geometry  of  7AZI  employing  the  AMI 
Hamiltonian. 

2- color  TOFMS  of  a)  7AZI(CH4)i  and  (b)  7AZI(Ar)i.  Wavelength  for  the  ionizing 
laser  is  3200  A. 

Optimum  geometries  for  7AZI(Ar)i  and  7AZI(CH4)i  obtained  with  a  potential  energy 
calculation.  Both  PM3  and  AMI  7AZI  structures  and  charges  result  in  identical  cluster 
geometries.  Also  shown  are  cluster  binding  energies. 

1 -color  TOFMS  of  a)  7AZI(H20)i  and  (b)  7AZI(H20)2-  The  peaks  between  33200  cm* 
1  and  33400  cm'1  are  2-3  times  larger  than  shown. 

MOPAC  5  results  for  the  local  hydrogen  bonding  structure  of  7AZI(H20)i  using  (a) 
PM3  and  (b)  AMI  Hamiltonians. 

The  potential  energy  calculation  results  for  7AZI(H20)i  using  (a)  PM3  molecular 
structures  and  charges  and  (b)  AMI  structures  and  charges.  The  percentages  listed  in  the 
Figure  refer  to  the  number  of  times  each  structure  is  found  starting  from  random  initial 
molecular  positions  and  orientations.  Also  shown  are  the  cluster  binding  energies  and  the 
percentages  of  the  outcomes. 

1 -color  TOFMS  of  7AZI  clustered  with  D2O.  Detection  mass  channels  are  indicated  in 
each  spectrum. 

1 - color  TOFMS  of  7AZI  clustered  with  D2O,  detected  at  mass  channel  (118+37).  (See 
text  for  a  discussion  of  the  doublet  (A,  C)  and  singlet  features.) 

2- color  TOFMS  of  a)  7AZI(CH30H)i  and  b)  7AZI(C2H50H)i.  Wavelength  of  the 
ionizing  laser  is  3225  A. 

The  potential  energy  calculation  result  for  the  geometry  of  7AZI(CH30H)i.  PM3 
molecules  are  used. 

2-color  TOFMS  of  a)  7AZI(NH3)i  and  b)  7AZI(NH3)2-  Wavelength  of  the  ionizing 
laser  is  3244  A. 


21 


Figure  15:  The  potential  energy  calculation  results  (based  on  PM3  molecules)  for  (a)  7AZI(NH3)i 
and  (b)  7AZI(NH3)2. 

Figure  16:  Comparison  of  the  vdW  modes  built  on  the  origins  of  various  clusters,  (a)  7AZI,  (b) 
7AZI(CH30H)i,  (c)  7AZI(Ar)i,  (d)  7AZI(CH4)l,  (e)  7AZI(H20)i,  (0  7AZI(H20)2, 
(g)  7AZI(H20)3,  (h)  7AZI(NH3)1,  0)  7AZI(NH3)2- 
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